Post-Newtonian accurate parametric solution to the dynamics of spinning compact 
binaries in eccentric orbits: The leading order spin-orbit interaction 
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We derive Keplerian-type parametrization for the solution of post-Newtonian (PN) accurate 
conservative dynamics of spinning compact binaries moving in eccentric orbits. The PN accurate 
dynamics that we consider consists of the third post-Newtonian accurate conservative orbital dy- 
namics influenced by the leading order spin effects, namely the leading order spin-orbit interactions. 
The orbital elements of the representation are explicitly given in terms of the conserved orbital en- 
ergy, angular momentum and a quantity that characterizes the leading order spin-orbit interactions 
in Arnowitt, Deser, and Misner-type coordinates. Our parametric solution is applicable in the fol- 
lowing two distinct cases: (i) the binary consists of equal mass compact objects, having two arbitrary 
spins, and (ii) the binary consists of compact objects of arbitrary mass, where only one of them is 
spinning with an arbitrary spin. As an application of our parametrization, we present gravitational 
wave polarizations, whose amplitudes are restricted to the leading quadrupolar order, suitable to 
describe gravitational radiation from spinning compact binaries moving in eccentric orbits. The 
present parametrization will be required to construct 'ready to use' reference templates for gravita- 
tional waves from spinning compact binaries in inspiralling eccentric orbits. Our parametric solution 
for the post-Newtonian accurate conservative dynamics of spinning compact binaries clearly indi- 
cates, for the cases considered, the absence of chaos in these systems. Finally, we note that our 
parametrization provides the first step in deriving a fully second post-Newtonian accurate 'timing 
formula', that may be useful for the radio observations of relativistic binary pulsars like J0737-3039. 
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I. INTRODUCTION 

Inspiralling compact binaries of arbitrary mass ratio, 
consisting of black holes and neutron stars, are the most 
plausible sources of gravitational radiation for the com- 
missioned ground based and proposed space based laser 
interferometric gravitational-wave detectors [1]. The 
possibility of detecting gravitational radiation from these 
sources are based on the following astrophysical and data 
analysis considerations. The detailed astrophysical in- 
vestigations indicate that the inspiral of compact bina- 
ries should be observed by gravitational-wave detectors 
[2, 3]. The temporally evolving gravitational wave polar- 
izations, h + (t) and h x (t), associated with the inspiralling 
compact binaries in quasi-circular orbits are known very 
precisely [4, 5]. This allows data analysts to employ 
the optimal method of matched filtering to extract these 
weak gravitational wave signals from the noisy interfero- 
metric data [6]. 

It should be noted that the construction of the 'ready 
to use search templates', namely above mentioned h+(t) 
and h x (t), is possible mainly due to the fact that the 
orbital dynamics, during the binary inspiral, is well de- 
scribed by the post-Newtonian (PN) approximation to 
general relativity. The PN approximation to general rel- 
ativity allows one to express the equations of motion for 
a compact binary as corrections to Newtonian equations 



'Electronic address: C.Koenigsdoerller@uni-jena.de 
t Electronic address: A.Gopakumar@uni-jena.de 



of motion in powers of (v/c) 2 ~ GM/(c 2 R), where v, M, 
and R are the characteristic orbital velocity, the total 
mass and the typical orbital separation of the binary. At 
present, the orbital dynamics of non-spinning compact 
binaries are explicitly known to 3.5PN order, which gives 
(v/c) 7 corrections to Newtonian equations of motion. 
These lengthy high PN corrections were obtained, for the 
first time, in Arnowitt, Deser, and Misner (ADM)-type 
coordinates [7-10]. Independently, they were later com- 
puted in harmonic coordinates and, as expected, found 
to be in perfect agreement [11-16]. Further, we point 
out that the very recent determination of the highly de- 
sirable 3.5PN (relative) accurate phase evolution for the 
gravitational wave polarizations was possible, mainly be- 
cause of employing techniques used in the computations 
of 3.5PN accurate equations of motion for the dynamics 
of compact binaries [9] . 

However, the effect of spins and orbital eccentricity on 
the compact binary dynamics and the associated h+(t) 
and hx{t) are not so accurately determined. The dy- 
namics of spinning compact binaries are determined not 
only by the orbital equations of motion for these objects, 
but also by the precession equations for the spin vectors 
themselves [17]. In fact, for the dynamics of spinning 
compact binaries only the leading order contributions to 
spin-orbit and spin-spin interactions are well understood 
[18, 19]. It is therefore quite customary, in the existing 
literature, to consider only the leading order contribu- 
tions to spin-orbit and spin-spin interactions, while ex- 
ploring the dynamics of spinning compact binaries [20] 
or the influence of spins on the phase evolution and am- 
plitude modulations of the gravitational waveforms [21]. 
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However, note that the equations of motion for spinning 
compact binaries, which include first post-Newtonian 
corrections to the leading order spin-orbit contributions 
are also available [22]. The leading order spins effects on 
the gravitational wave polarizations, mainly for compact 
binaries in quasi-circular orbits, were also explored in 
great detail [23] . As far as the orbital eccentricity is con- 
sidered, the temporally evolving h+(t) and h x (t), whose 
amplitudes are Newtonian accurate and the orbital evo- 
lution is exactly 2.5PN accurate, was recently obtained 
[24]. However, the inputs required to compute the cru- 
cial secular phase evolution for compact binaries in inspi- 
ralling eccentric orbits are also available to 2PN order be- 
yond the leading quadrupole contributions [24, 25]. Fur- 
ther, very recently, a Keplerian-type parametrization for 
the solution of 3PN accurate equations of motion for two 
non-spinning compact objects moving in eccentric orbits 
was obtained [26]. It is therefore desirable to obtain a 
Keplerian-type parametric solution to the dynamics of 
compact binaries, when spin effects are not neglected. 

In this paper, we obtain a Keplerian-type parametric 
solution to the PN accurate dynamics of spinning com- 
pact binaries moving in eccentric orbits, when the dy- 
namics includes the leading spin-orbit interaction. Our 
parametric solution is applicable in the following two dis- 
tinct cases: (i) the binary is composed of two compact 
objects of equal mass and two arbitrary spins, and (ii) 
the two objects have unequal masses but only one object 
has a spin. The present parametrization can easily be 
added to the recently obtained 3PN accurate generalized 
quasi-Keplerian parametrization for the solution of the 
3PN accurate equations of motion for two non-spinning 
compact objects moving in eccentric orbits [26]. This will 
provide, for the first time, almost analytic, 3PN accurate 
parametric solution to the dynamics of spinning com- 
pact binaries moving in eccentric orbits, when the spin 
effects are restricted, for the time being, to the leading 
order spin-orbit contributions [27]. As an application of 
the parametrization, we present explicit expressions for 
the spin-orbit modulated h + {t) and h x (t), whose ampli- 
tudes are restricted, for convenience, to the Newtonian 
order and the orbital motion is restricted to the above 
mentioned conservative PN order. Our results generalize 
and improve, for the two cases considered, a restricted 
parametrization for the orbital motion of compact bina- 
ries with leading order spin-orbit coupling [28]. 

Our parametrization will have the following possible 
applications. The construction of 'ready to use' search 
templates, namely h+(t) and h x (t) 7 for compact binaries 
in inspiralling orbits requires a detailed knowledge of the 
PN accurate orbital dynamics and a PN accurate descrip- 
tion for hf?, the transverse-traceless (TT) part of the ra- 
diation field. As mentioned earlier, time evolving h + (t) 
and h x (t) for non-spinning compact binaries of arbitrary 
masses moving in inspiralling eccentric orbits were ob- 
tained in Ref. [24]. The parametrization presented here, 
along with [26], will be required to construct PN accurate 
'ready to use' reference templates for spinning compact 



binaries moving even in inspiralling eccentric orbits. 

There is an ongoing debate about whether or not the 
spinning compact binaries, of arbitrary mass ratio, mov- 
ing in circular or eccentric orbits under PN accurate dy- 
namics will exhibit chaotic behavior [20]. The question 
of chaos was originally motivated by the observation that 
the equations describing the motion of spinning test par- 
ticle in Schwarzschild spacetime allow chaotic solutions 
[29]. The issue is usually addressed by solving numer- 
ically PN accurate equations describing the dynamics 
of spinning compact binaries to compute certain gauge 
dependent quantities like fractal basin boundaries and 
Lyapunov exponents. The semi-analytic parametrization 
presented here, along with the future extensions, should 
be very effective in analyzing the PN accurate dynam- 
ics of spinning compact binaries and to explore if astro- 
physically realistic binary inspiral will be chaotic or not. 
Further, the parametrization for the cases considered im- 
plies that the associated PN accurate conservative binary 
dynamics can not exhibit chaos. 

Finally, we note that our parametrization may be use- 
ful to analyze the high precision radio-wave observa- 
tions of relativistic binary pulsars like J0737-3039 [30]. 
These radio-wave observations employ a 'timing for- 
mula', which gives the arrival times at the barycenter of 
the solar system for the electromagnetic pulses emitted 
by a binary pulsar. The 'timing formula' is important 
for obtaining astrophysical informations from the com- 
pact binary as well as to test general relativity in strong 
field regimes [31, 32]. The heavily employed 'timing for- 
mula', usually incorporates 1PN accurate orbital motion, 
leading order secular effects due to gravitational radia- 
tion reaction, spin-orbit coupling as well as 2PN (secu- 
lar) corrections to the advance of periastron [31, 33-35]. 
The parametric solution obtained here will be required to 
obtain a 2PN accurate 'timing formula', where all spin- 
orbit and 2PN corrections are fully included. This im- 
plies a 'timing formula', which includes not only 2PN 
order secular effects, but all 2PN order periodic terms, 
which may leave some potentially observable signature 
[28, 36, 37]. 

The paper is structured in the following manner. In 
Sec. II, we exhibit and explain, in detail, the total Hamil- 
tonian describing the conservative dynamics, whose para- 
metric solution we are seeking in this paper. In Sec. Ill, 
we summarize the Keplerian-type parametrization that 
describes the solution to the 3PN accurate equations of 
motion for compact binaries, in eccentric orbits, when 
spin effects are ignored. Sec. IV describes the determina- 
tion of a Keplerian-type parametrization for the conser- 
vative dynamics of spinning compact binaries, when spin 
effects are restricted to the leading order spin-orbit in- 
teraction. The parametric solutions, presented in Sec. Ill 
and Sec. IV, are combined in Sec. V, to obtain PN ac- 
curate Keplerian-type parametrization for the conserva- 
tive dynamics of spinning compact binaries moving in 
eccentric orbits. As an application to our PN accurate 
parametric solution, we present explicit analytic expres- 
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sions for the gravitational wave polarizations, h+(t) and 
h x (t), for spinning compact binaries in eccentric orbits 
in Sec. VI. The temporal evolution of above h + (t) and 
h x (t), whose amplitudes are Newtonian accurate, will be 
governed, for the time being, by the conservative PN dy- 
namics discussed in this paper. Finally, in Sec. VII, we 
summarize our results and discuss its future extensions 
and applications. 



II. DYNAMICS OF COMPACT BINARIES 
WITH LEADING ORDER SPIN ORBIT 
INTERACTION 

In this paper, as mentioned earlier, we are interested 
in the 3PN accurate conservative dynamics of spinning 
compact binaries, when spin effects are, for the time be- 
ing, restricted to the leading order spin-orbit interaction. 
The dynamics is fully specified by a post-Newtonian ac- 
curate Hamiltonian H, which may be symbolically writ- 
ten as 

H = 7Yn + Wipn + W2PN + ^3PN + Hso i (2.1) 

where 7Yn, T^ipn, W2PN and Tispn respectively are the 
Newtonian, first, second and third post-Newtonian con- 
tributions to the conservative dynamics of compact bi- 



naries, when the spin effects are neglected. The leading 
order spin-orbit coupling to the binary dynamics is given 
by Ti-so . We determine the parametric solution to the dy- 
namics, given by Eq. (2.1), in the following way. First, we 
obtain a parametric solution to the PN accurate conser- 
vative dynamics, neglecting the effects due to spin-orbit 
coupling. We then compute a parametric solution to the 
conservative dynamics, given by 7^n + ^so- The second 
step is consistent in a post-Newtonian framework as the 
spin-orbit interactions are restricted to the leading or- 
der. In the final step, we consistently combine these two 
parametrizations, to obtain a PN accurate parametric so- 
lution to the conservative compact binary dynamics, as 
specified by TL. 

In the present work, we employ the following 3PN ac- 
curate conservative Hamiltonian, in Arnowitt, Deser, and 
Misner (ADM)-type coordinates [38, 39], supplemented 
with a contribution describing the leading order spin- 
orbit coupling. We work, as is usual in the literature 
[26], with the following 3PN accurate conservative re- 
duced Hamiltonian, H = 7i//i, where the reduced mass 
of the binary is given by jx = mim2/M, mi and m 2 
being the individual masses and M = mi + m 2 is the 
total mass. The 3PN accurate (reduced) Hamiltonian, in 
ADM-type coordinates and in the center-of-mass frame, 
with the leading order spin-orbit contributions reads 



H(r,p,Si,S 2 ) = H N (r,p) + Hi PN (r,p) + H 2PN (r,p) + H 3PN (r,p) + H so (r,p,Si,S 2 ) , (2.2) 
where the Newtonian, post-Newtonian and spin-orbit contributions are given by 



H N (r,p) 



p 2 1 



H 1P N(r >P ) = 1 far, 1) (p 2 ) 2 - \ [(3 + v)p 2 + »(» ■ P?} \ + 



H 2PN (r,p) = ^ (1 - 57] + W) (p 2 



(5 - 20r) - 3?? 2 ) (p 2 ) 2 - 2r/ 2 (n • p) 2 p 2 - 3t7 2 (ti • p) L 



1 [(5 + 8 V )p 2 + 377(71 ■ p) 2 ] \ - j(l + 3r ? ) 1 



H 3 PN(r,p) 



2 

1/1 
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r 

4 1 [ 
+ 16 



(2.3a) 

(2.3b) 

1 

r 

(2.3c) 
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16 



(-27 + 1367 7 + 109?7 2 ) (p 2 



) 2 + -L(17 + 3077)77(71 ■ pfp 2 + ^(5 + 43r7)r7(n • pf 



[y^[- 600 + ( 3n2 - 1340 ) r ' - 552? ?V - ll ( 340 + 3 ^ 2 + H2 77 ) 77(n ■ p) 2 | ^ 



-l[l2+ (872 -63^)77]^ 
1 



H so (r,p, Si,S 2 ) = -^(rxp) • S eB 



(2.3d) 
(2.3e) 
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In the above equations r = TZ/(GM), r — \r\, n = r/r 
and p = V/fi, where TZ and V are the relative separation 
vector and its conjugate momentum vector, respectively. 
The finite mass ratio r\ is given by 77 = /Lt/M . The spin- 
orbit coupling involves the effective spin S e g, given by 

S eS = 5 1 S 1 + 5 2 S 2 , (2.4) 

where 

*- 2 "( i+ £t)^ + ?( i+ ^)^ (2 - 5b » 

The reduced spin vectors Si and S 2 are related to the 
individual spins Si and S 2 by S\ = Si/ (nGM) and 
S 2 = S 2 /{jiGM), respectively. We recall that TZ, V, 
Si and S 2 are canonical variables, such that the or- 
bital variables commute with the spin variables, e.g. see 
Ref. [17, 19]. 

The above definition of the effective spin S e s is iden- 
tical to the quantity £, introduced in Ref. [28], and re- 
lated to SWItd, introduced in Ref. [19], by 2S , e ff|xD = 
GM 2 S eS . 

The PN corrections, associated with the motion 
of non-spinning compact binaries, are compiled using 
Refs. [9, 40]. The spin-orbit contribution to H, avail- 
able in Ref. [35], employs the spin-supplementary condi- 
tion (SSC) of Pryce, Newton and Wigner [41, 42]. The 
SSC, which defines the central world line of each spinning 
body, adopted in this paper is given by 

2S M + - c S l3 W = (ij = 1, 2, 3) , (2.6) 

where <S M „ is the spin-tensor and J7 M the 4-velocity of 
the center of mass of the body. The spin 4-vector S a is 
defined by 

S a = —ea^uU^ (a, p,fi,v = 0,1, 2, 3), (2.7) 

where Eap^ is the Levi-Civita tensor with £0123 = +1- 
In the rest frame of the body (£/ M = cSq ) 

S a = (0,S 23 ,S 31 ,S 12 ) = (0,5) (2.8) 

holds. The spins, «Si and 52, of the compact objects 
may be given in terms of their moments of inertia and 
angular velocities of proper rotations, I a and fl a (a = 
1,2) respectively, as 

51 =X 1 fl 1 , (2.9a) 

5 2 =l 2 n 2 . (2.9b) 

Using the above definitions, it is possible to determine at 
what PN order the leading order spin-orbit interaction 
will manifest. It is easy to see, using Eq. (2.3), that spin 



effects enter the dynamics formally at 1PN order. How- 
ever, for compact objects, the spin is roughly given by 
S ~ m co r co v spin ~ G m 2 v spm / c 2 , where m co ,r co and 
v spm are typical mass, size and rotational velocity of 
the spinning compact object. This implies, for moderate 
V spm values, that the leading order spin orbit coupling 
will enter the compact binary dynamics at 0(l/c 4 ), i.e. 
at 2PN order. It is interesting to note that if v spin ~ c, 
then the leading order spin-orbit interaction will numer- 
ically manifest at 1.5PN order. However, in that case the 
binary dynamics should be dominated by the currently 
neglected PN corrections to H^o- In this paper, we will, 
for the sake of clarity and simplicity, assume that the 
spin-orbit interaction influence the compact binary dy- 
namics at 2PN order. We would like to point out that the 
present derivation of the Keplerian-type parametrization 
for the PN accurate dynamics of spinning compact bina- 
ries does not really care if the spin-orbit coupling mani- 
fests at 1.5PN or 2PN order. 

In the next section, we will display and explain the 
recently obtained Keplerian-type parametric solution to 
the 3PN accurate conservative orbital motion of two com- 
pact objects, of negligible proper rotations, moving in 
eccentric orbits [26]. 



III. KEPLERIAN-TYPE PARAMETRIZATION 
FOR THE ORBITAL MOTION OF COMPACT 
BINARIES WITH NEGLIGIBLE PROPER 
ROTATION 

Using the fully determined 3PN accurate conserva- 
tive Hamiltonian, in ADM-type coordinates and in the 
center-of-mass frame, available in Refs. [9, 40], very re- 
cently a Keplerian-type parametrization for the orbital 
motion of compact binaries in eccentric orbits was de- 
rived [26]. The derivation of the above parametric so- 
lution, usually referred to as the 3PN accurate general- 
ized quasi-Keplerian parametrization, crucially depends 
on the following important points. First, the 3PN ac- 
curate conservative Hamiltonian is invariant under time 
translation and spatial rotations, implying the existence 
of 3PN accurate (reduced) energy E = H and (reduced) 
orbital angular momentum L = r x p, for the binary in 
the center-of-mass frame. In particular, the conservation 
of L indicates that the motion is restricted to a plane, 
namely the orbital plane, and we may introduce polar 
coordinates such that r = r(cos tp, sinyj). Finally, the 
Hamiltonian equations of motion governing the relative 
motion of the compact binary, namely the 3PN accurate 
expressions for f = dr/dt and <p — dip/dt — where t de- 
notes the coordinate time scaled by GM — in terms of 
E,L = \L\ and r, are polynomials of degree seven in 1/r. 

The third post-Newtonian accurate generalized quasi- 
Keplerian parametrization for compact binaries moving 
in eccentric orbits, in ADM-type coordinates and in the 
center-of-mass frame, is given by 
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r = a r (1 — e r cos u) , 



<f — <Po 



V c 4 



i / 1 i \ • . (9u , 9u\ , ^ , f fit , fet\ . . ? 6t . „ . ftet • 
1 — n (t — to) = u — e t sin u + I — j- H — ^- 1 (v — u) + I —r- H — g- I sin u H — g- sin 2i> H =- sin 3v , 



(1 + k ) v + ( i*f + % ) S m2v + ( ^ + ^ sin3« + ^ sin4 V + % sin5 
V c 4 c 6 / V c 4 c 6 / c 6 c 6 



where v = 2 arctan 



l + e ¥ 
1 — e u 



1/2 



tan ■ 



r 



(3.1a) 
(3.1b) 

(3.1c) 
(3.1d) 



The PN accurate orbital elements a r , el 



k, and 



e v , and the PN orbital functions g it , get, fit, fet, iet, 
het, f4 V , fetp, 34 V , ge v , ie v , and h Gtp expressible in terms 
of E, L and rj are obtainable from Ref. [26] and will be 
explicitly displayed in Sec. V, after including the effects 
due to the leading order spin-orbit interactions. 

Let us make a few comments about the parametriza- 
tion. The radial motion is uniquely parametrized by 
Eq. (3.1a) in terms of some PN accurate semi-major axis 
a r , radial eccentricity e r and the eccentric anomaly u. 
The angular motion, described by Eq. 3.1c, is specified 
by the true anomaly v, angular eccentricity e v , advance 
of the periastron k and some 2PN and 3PN order orbital 
functions. The explicit time dependence is provided by 
the 3PN accurate 'Kepler equation', namely Eq. 3.1b, 
which connects the mean anomaly I, and hence the co- 
ordinate time, to eccentric and true anomalies u and v, 
respectively. The PN accurate 'Kepler equation' also in- 
cludes some PN accurate mean motion n, time eccentric- 
ity et and some 2PN and 3PN order orbital functions. It 
should be noted that all these eccentricities, e r , e v and et 
are connected to each other by PN accurate expressions 
in terms of E, L and r], given by Eq. (21) in Ref. [26]. This 
implies that it is possible to specify, as in the Newtonian 
case, a PN accurate eccentric orbit in terms of, for exam- 
ple, a r and e r . Finally, we recall that the Keplerian-type 
parametrization for compact binaries in eccentric orbits 
at first post-Newtonian order was obtained in Ref. [43]. 
Rs extension to 2PN order was derived in Refs. [35, 44]. 

In the next section, we will improve the above 
parametrization by including the effects due to the lead- 
ing order spin-orbit interactions. 



IV. QUASI-KEPLERIAN PARAMETRIZATION 
OF THE CONSERVATIVE DYNAMICS OF A 
SPINNING COMPACT BINARY — 
NEWTONIAN MOTION AUGMENTED BY THE 
FIRST-ORDER SPIN-ORBIT INTERACTION 

In this section, we incorporate in to the parametriza- 
tion, obtained in the previous section, the effects due to 
the leading order spin-orbit coupling. The straightfor- 
ward way of implementing that is to analyze the orbital 
dynamics, which is simply the Newtonian one augmented 



by the leading order spin-orbit interaction. Let us first 
describe, in the Hamiltonian framework, the binary dy- 
namics we are interested in. 



A. The Newtonian binary dynamics with leading 
order spin effects 

We investigate the binary dynamics, given by the re- 
duced Hamiltonian -Hnso , where -Hnso = + H^o ■ 



H P 1 

-Hnso — -7, 

Z r 



1 



LS 



eff • 



.1) 



The above Hamiltonian prescribes, via the Poisson brack- 
ets, the following evolution equations for the reduced an- 
gular momentum vector L = r x p, and the individual 
reduced spin vectors Si and 



^ = {L,H NBO } = 



1 



S e ff x L . 



Si 
C 5 2 



L x Si 



L x S 2 



(4.2a) 
(4.2b) 
(4.2c) 



where {...,...} denotes the Poisson brackets. Note that 
dL/dt, given by Eq. (4.2a), gives only the precessional 
motion of the orbital plane. We still need to obtain the 
equations describing the orbital motion. Before we go 
onto derive those equations, let us consider the conserved 
quantities associated with -Hnso- 

The reduced energy E = Hnso is conserved because 
of 9 t 7?NSO = 0. Though L is not conserved, it is fairly 
easy to show that its magnitude L = |L| is conserved: 



:i'(S«ffXi) = 0. (4.3) 



dL 2 _ d 
~~dT ~ di 



Similarly, below we show that the magnitudes of the spins 
are also conserved: 



dS 2 
dt 

,'2 



dt 



(Si • Si) 



ITT "Si ■ (L x S a 



26 2 



, (4.4a) 
S 2 -(ixS 2 )=0. (4.4b) 



Note that Eqs. (4.2) indicate that S — —L, where S = 
Si + S 2 is the total spin reduced by fj,GM, the reduced 
total spin. We now define the conserved (reduced) total 
angular momentum vector 



J 



(4.5) 



The above expression times fiGM gives the total angular 
momentum namely = 1Z x "P + S. Moreover, 
J is conserved both in its magnitude and direction as 
dJ/dt = and dJ/dt = 0, where J = \ J\. 

The fact that dJ/dt = allows us to introduce a ref- 
erence orthonormal triad (ex,ey,ez), such that ez is 
along the fixed vector J. This gives 



J = Je 



z ■ 



(4.6) 



Further, for later use, we choose the line-of-sight unit 
vector N to lie in the ey-e^-plane (see Fig. 1). This is 
always possible and reasonable because of the degree of 
freedom associated with the choice of ex and ey ■ These 
unit vectors ex and ey span the invariable plane which 
is the plane perpendicular to J. 

We now introduce highly advantageous spherical (po- 
lar) coordinates (r, 6, <j>) and the associated orthonormal 
triad (n,eg,e,p), where 9 is the angle between ez and 
n and 4> is the azimuthal angle defined in the invariable 
e^-ey-plane (see Fig. 1), such that 



sin 6 cos <p ex + sin u sin <p ey + cos uez 



(4.7) 



The relative separation vector r and its time derivative 
r = dr/dr, written in the orthonormal triad (n,eg,e^,) 
reads 



r = rn , 

r = rn +r6 e$ + r sin 9<f> e^ 



(4.8a) 
(4.8b) 



In spherical coordinates, the linear momentum p and its 
magnitude may be expressed as 



p = p r n + peee + p^e^ , 

2 j 
- 

L 2 



p 2 =pl+p 2 e +pl = {n- p) 2 + (n x p) 2 



Pr + ~n ■ 



(4.9a) 



(4.9b) 



With the help of these relations and E = i^NSO, the 
components of the linear momentum p can be expressed 
using the quantities E, L and (L ■ S e g) as 



p 2 r =2E- 



2 L 2 2{L ■ S eS ) 



P<t> 



L z 

r sin( 

-2 



2 ^ 2 1 / j- 2 ^ Z 

Pd = ~ Pd> = ~ I -k . 2 n 

r z v r z \ sin 6 



(4.10a) 
(4.10b) 
(4.10c) 



where L z = e z ■ L. Note that in general L z is not 
conserved. 



The components of r in the spherical polar coordi- 
nates, which also define the orbital equations of motion 
associated with ^nso, read 



r = n ■ r = p r , (4.11a) 
4 = e e ■ r = V e + , (4.11b) 

(4.11c) 



eg ■ S 



eff 



In the above equations, we have used Hamilton's equation 
r = dH^so/dp, n x eg = e^ and nxe^ = —eg. 

Summarizing, we note that the dynamics of the bi- 
nary system, described by the Hamiltonian -Hnso, is also 
uniquely determined by the evolution equations given by 
Eqs. (4.2) and (4.11). Before trying to find parametric 
solutions to these equations, let us point out several fea- 
tures of the above binary dynamics. 

It is easy to write down the precessional equations for 
the (reduced) total spin S and the effective spin S e g as 



S = Si + S 2 



L x S e ff = L . 
1 



(4.12a) 



S eff = 5 1 S 1 + S 2 S 2 = x {5(Si + 5 2 S 2 ) . (4.12b) 

It is clear that the motions of S and S e s are more com- 
plicated than those for S\ and S 2 . We deduce that the 
magnitudes of S and S^ff are not preserved in general 
and satisfy following equations 



dS 2 
dt 



-^S.( S ,x S2 ), 



(4.13a) 



djk-_ h^m^± rn L . {SixS2h (413b) 



dt 



Ac 2 r 3 



We note that the quantity (L-S e g) is exactly preserved 
since 

|(r-5-)-f-»- + X.^-0. (414, 

Finally, let us consider the evolution of Lz. Using 
L = Jez — S on the r. h. s. of Eq. (4.2a) we obtain 



dL z 3VT^4^ 



dt 



2c 2 ? 



e z -(S 1 xS 2 ). 



(4.15) 



The above equations imply that, in general, Lz is not 
conserved. However Lz along with S and 5 e ff are con- 
served for the following two cases. In the first instance, 
case (i), we require the binary to have equal masses, but 
with arbitrary spins Si and S 2 (mi = m 2 Tj = 1/4 
^1 = ^2). In the second instance, case (ii), the binary may 
have arbitrary masses, but only one of them is spinning 
(mi 7^ m 2 , Si ^ or S 2 ^ 0). 

In this paper, when we solve the differential equations, 
Eqs. (4.2) and (4.11), we restrict ourselves, for the time 
being, to the above two cases. In these cases, the effective 
spin iSeff and the reduced total spin S are related by 



>eff 



S1S1 + S 2 S 2 = XsoS , 



(4.16) 
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where the newly defined mass dependent coupling con- 
stant Xso is given by 



Xso • 



5\ = 82 = 7/8 for (i), the equal-mass-case, 



Si or 82 



for (ii), the single-spin-case. 



We would like to emphasize that in the case (i) the re- 
duced total spin S is the sum of two arbitrary spins Si 
and S2. However, in the case (ii) S is either Si or S^. 

These observations allow us to present the precessional 
equations for L and S in the following compact form for 
the special cases as 



dL 

IE 

dS 

dt 



(4.17a) 
(4.17b) 



The Eqs. (4.17) imply that the (reduced) orbital angular 
momentum L and the (reduced) total spin S, for the 
two cases considered, precess about the fixed vector J = 
Jez at the same rate with an (instantaneous) precession 
frequency given by 



Xs 



(4.18) 

u— / " 

Note that uj p can only be defined in a physically mean- 
ingful way when spin effects are included in the binary 
dynamics. In the non-spinning limit, the cross prod- 
ucts appearing on the r. h. s. of Eqs. (4.17) vanish, 
leaving the (reduced) orbital angular momentum L as a 
constant vector. It is interesting to observe that in the 
non-spinning limit, neither lu p nor x so go to zero. 

In the next subsection, we will find a parametric so- 
lution to the dynamics of spinning compact binaries in 
eccentric orbits, given by Eqs. (4.11) and (4.17). 



B. The Keplerian type parametrization associated 
with the Hamiltonian -Hnso 

We start by considering the radial motion, governed 
by Eq. (4.11a), which reads 



f 2 = 2E 



L 2 2(L ■ S eB ) 



(4.19) 



We obtain parametric solution to the above equation 
by following exactly the same procedure described in de- 
tail in Sec. Ill of Ref. [26]. The radial motion is described 
by 



r = a r (l — e r cosw) , 
I = n(t — to) = u — e* sin u , 



(4.20a) 
(4.20b) 



where u and I are the eccentric and mean anomalies. The 
orbital elements, explicitly given in terms of E, L and 
L ■ S e ff, are 



1 /, r,L &ef[ E 
CLr = 1 — 2 - - 

2E \ L 2 c 2 



(4.21a) 



e 2 = 1 + 2EL 2 



n = (-2Ef' 2 



1 + 2EU 



8(1 + EI? 



L ■ Sgg E 

L 2 c 2 



I L ■ S e ff E 
L 2 c 2 



(4.21b) 
(4.21c) 
(4.21d) 



Eq. (4.20b) gives the 'Kepler equation', modified by 
the spin-orbit interaction, which connects the eccentric 
anomaly u to the coordinate time. The crucial require- 
ments to determine the above parametrization are the 
conservation of E, L and L ■ S e ff. This allows us to treat 
(dr/dt) 2 as a polynomial in 1/r with constant coefficients. 

Let us now find parametric solution to the angular 
parts of the orbital equations of motion, Eqs. (4.11) com- 
bined with Eq. (4.16), written in the form 



rB 



Ys 



Pe [I 



Xso J sm9 



(4.22a) 
(4.22b) 



For the cases of interest, where Lz is a constant, we may 
introduce a constant angle between ez and L such that 
L z = LcosG (see Fig. 1) and Eqs. (4.10b) and (4.10c) 
simplify to 



L cos O 



r smt 



cos 2 
sin 2 6 



(4.23a) 
(4.23b) 



The efficient way of determining solution to the 
Eqs. (4.22) begins by expressing the radial vector r in 
terms of an orbital triad. This is achieved by connecting 
via two rotations the reference triad (ex,ey,ez) to an 
orbital triad defined by k). In our orbital orthonor- 
mal triad, the unit vector fe, given by k = L/L, and it 
along with ez defines 



% = 



ez x k 
\e z x fe| 



(4.24) 



Note that \ez X fe| = sin©, which is indeed always posi- 
tive because of < 9 < n. We note that the unit vector 
i gives the line of nodes, associated with the intersection 
of the orbital plane with the invariable plane (ex,ey). 
The corresponding inclination angle is 9. Since the above 
mentioned line of nodes vanishes in the non-spinning 
limit, the unit vector i is only defined when spin effects 
are included in the binary dynamics. 

The two rotations arise from the observation that the 
(instantaneous) position and orientation of the orbital 
plane with respect to the reference triad (ex, ey,ez) are 
defined by two angles: the longitude of the line of nodes 
T, (0 < T < 2tt) and the inclination of the orbital plane 
with respect to the invariable ex-ey-plane 9, (0 < < 
7r). In terms of rotational matrices, we have 



TO 
cos 9 sin 9 
. — sin 9 cos 9 , 



cos T sin T 0\ (ex 
— sin T cos T ey 
1/ \e z . 
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In the new orbital orthonormal triad (i,j,k), the rel- 
ative separation vector r is given by 



r = r cos tpi + r sin ip j , 



(4.25) 



where the unit vectors (i,j,k) in terms of the reference 
triad (ex,ey,ez) are explicitly given by 

i = cos T ex + sin T ey , (4.26a) 

j = — cos 9 sin T ex + cos 9 cos T ey + sin 9 ez , 

(4.26b) 

k = sin 9 sin T ex — sin 9 cos T ey + cos 9 ez , (4.26c) 

The geometrical interpretations of the newly introduced 
angles and basic vectors are illustrated in Fig. 1. 



k = L/L 



Line of sight 




Figure 1: The binary geometry and interpretations of vari- 
ous angles appearing in this section. Our reference frame is 
(ex,ey,ez), where the basic vector ez is aligned with the 
fixed total angular momentum vector J = L + S , such that 
J = Jez. The invariable plane (ex, ey) is perpendicular 
to J. Important for the observation is the line-of-sight unit 
vector N from the observer to the source (compact binary). 
We may have, by a clever choice of ex and ey, the line-of- 
sight unit vector N in the e.\ — ez-plane. k = L/L is the 
unit vector in the direction of the orbital angular momentum 
L, which is perpendicular to the orbital plane. The constant 
inclination of the orbital plane with respect to the invariable 
plane is 8, which is also the precession cone angle of L around 
J. The orbital plane intersects the invariable plane at the line 
of nodes i, with the longitude T measured in the invariable 
plane from ex. T is also the phase of the orbital plane preces- 
sion. The orbital plane is spanned by the basic vectors 
where j = k x i. 



The comparison of r, given by Eqs. (4.8a) and (4.7), 
with the one in the new angular variables, Eq. (4.25) with 
Eqs. (4.26), implies the following transformations for the 
angular variables: 

!cos 9 = sin ip sin 9 
sin(0 — T) sin 9 = sin ip cos 9 
cos(</> — T) sin 9 = cos p 

(4.27) 

A clever combination of the above angular transforma- 
tion equations and their corresponding time derivatives 
leads to the following relations for the angular velocities: 



= t 



cos 2 9 
sin 2 9 
cos 9 
sin 2 9 1/5 



<P 



(4.28a) 
(4.28b) 



where the over dot again means the derivative with re- 
spect to time t and 9 is treated as a constant angle. 

Using Eq. (4.28a) for 9 2 in Eq. (4.22a) with (4.23b) 
leads to 



L (, Xso 
<P = — I 1 ~ T" 



(4.29) 



Note that the above derivation requires 9 ^ 0, which is 
always satisfied for S ^ 0, and hence Eq. (4.29) is not 
defined when 9 = 0. We again invoke the procedure 
employed in Sec. Ill of Ref. [26] to obtain the following 
parametric solution for p: 



y- <A) = (l + k)v 

v = 2 arctan 



\ 1/2 

tan^ 
e v J 2 



(4.30a) 
(4.30b) 



where tpo is the value of <p at the time to. The quantity 
k is a measure of the advance of the periastron and e v is 



the so called angular eccentricity, 
in terms of E, L and L ■ S e g as 



They are expressible 



1 



c 2 i 2 



Xso - 3 



L ■ 5 e fr 
L 2 



1 + 2EL 2 -4 Xs o(l + 2 J Bi 2 )4 



4(3 + AEL' 



L ■ S e ff E 
L 2 c 2 



(4.31a) 



(4.31b) 



The parametrization for r, I and ip can also be ob- 
tained using the conchoidal transformations employed in 
Ref. [33] and found to be in total agreement. 

We now move on to derive the time evolution of T, the 
longitude of the line of intersection (denoted by the line 
of nodes unit vector i in Fig. 1), in a parametric way. 
Using Eq. (4.28b) for cf> with Eq. (4.29) in Eq. (4.22b) 
with Eq. (4.23a) we get the following differential equation 
for T: 



dT 



XsoJ 



(4.32) 
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which is also restricted to the spinning case (0 7^ 0), 
since we made use of Eq. (4.29) for ip. 

Note that the r. h. s. of the above equation is identical 
to u>p, given by Eq. (4.18). This is not surprising as 
the frequencies of the precessional motion for k = L/L 
and i should be identical due to Eqs. (4.24) and (4.26). 
Further, a close inspection of Fig. 1 reveals that the phase 
of the projection of k on to the invariable plane (ex, ey), 
as measured from ex, is given by T + 270°. 

To solve the differential equation (4.32) analytically, 
we divide it by dip/dt, given by Eq. (4.29), and deduce 
at the leading order 



dT _ XsqJ 1 
dip 



L err 



(4.33) 



The radial separation r appearing on the r. h. s. of 
above equation may be replaced by following expression 
in terms of v as shown below 



r = a(l — ecosu) 



L 2 



a(l-e 2 ) 



1 + e cos v 1 + e cos v 



(4.34) 



where a and e may be treated as the Newtonian accurate 
expressions for a r and e r . 

This allows us to write Eq. (4.33) — within the accu- 
racy needed — as 



dT XsoJ 



dv 



J L 3 



(1 + e cos v) 



(4.35) 



where we also used the fact that <p = v at the leading 
order. The Eq. (4.35) is easily integrated to obtain the 
parametric solution to T as 

T — To = . (v + e sin v) , (4.36) 

where we put v(to) = 0. 

There exists an alternate way to obtain dT/dt, given 
by Eq. (4.32). Noting that the time evolution for the unit 
vector k, defined by Eq. (4.26c), is solely given by that 
for T, we compute the time derivative of fe and it reads 



k = TsinO(cos Tex + sin Tey) 



(4.37) 



Using the fact that L = Lk and L = Lk in Eq. (4.17a), 
we arrive (again) at dT/dt, as given by Eq. (4.32). There- 
fore, the parametric solution for T readily defines a simi- 
lar solution for L = Lk, where k is given by Eqs. (4.26c) 
with (4.36). The parametric solution to L reads 

L = L(s'm sin T ex — sin cos T ey + cos ez) ■ 

(4.38) 

Finally, let us derive a parametric solution to the 
precessional motion of S, given by Eq. (4.17b). This 
parametric solution follows immediately by noting that 
S = J — L, where J = Jez and L = Lk. Using 
Eq. (4.26c) for fc, we deduce that 

S = L sin 0(— sin T)ex + L sin cos Tey 



+ ( J — LcosQ)ez , 



(4.39) 



The above equation, along with the parametric solution 
for T, describes, in a parametric way, the time evolution 
of S. 

We finally collect all the relevant equations, namely 
Eqs. (4.20), (4.30), and (4.36), and display below our 
parametric solution for the binary dynamics given by 



r = r cos ip % + r sin ip j , 

L = Lk, 

S = Jez — Lk , 



(4.40) 
(4.41) 
(4.42) 



where the basic vectors (i,j,k) are explicitly given by 

i = cos Tex + sin Tey , (4.43a) 
j = — cos sin Tex + cos cos Tey + sin ez , 

(4.43b) 

k = sin sin Tex — sin cos Tey + cos ez , (4.43c) 

The time evolution for the radial and angular variables 
is given by 



r = a r ( 1 — e r cos u) , 
n(t — to) = u — et sin u , 
<P - po = (1 + k)v , 

XsoJ , 



T-Tr 



2 L 3 



(v + esinw) , 



v = 2 arctan 



1 + e^V^ 2 , u 
tan — 



where the orbital elements are given by 



S E 

— 2x so cosa- — 



1 

2E 



S E 

1 + 2EL 2 + 8(1 + EL 2 ) Xso cos a--j 

L c z 



n = (-2E) 



3/2 



1 + 2EL 2 + 4y so cos a- , 



S E 

'Lc 

C 2 L 2 \X so ~ 3 Xs cosa — 



1 + 2EL 2 -A(l + 2EL 2 ) Xso ^ 



4(3 + AEL 2 )x S o cos a 



S E 



Lc 



2 ' 



(4.44a) 
(4.44b) 
(4.44c) 

(4.44d) 
(4.44e) 



(4.45a) 

(4.45b) 
(4.45c) 
(4.45d) 

(4.45e) 



(4.45f) 



As noted earlier, e is given by the Newtonian contribu- 
tion to e r . In Eqs. (4.45), we have replaced L ■ S e s by 
XsoLS cos a, where a is the constant angle between L 
and S and x so is given by 

J Si = 82 = 7/8 for (i), the equal-mass-case, 



Si or 5 2 



for (ii), the single-spin-case. 
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Further, we note that the constant angle is not a free 
variable and is given by one of the following relations: 



sin 9 



cos = 



Ssina 

L + S cos a 
J 



(4.46a) 
(4.46b) 



where the magnitude of the (reduced) total angular mo- 
mentum is given by J = (L 2 + S 2 + 2LS cos a) 1 / 2 . 

We are aware that the above parametrization does not 
lead to the classic Keplerian parametrization simply by 
putting S = 0. There are several arguments for this ap- 
parent lack of a simple non-spinning limit. Notice that 
T and <p are defined with respect to the line of inter- 
section, characterized by the line of nodes unit vector i. 
When the spin effects are neglected, the above line of 
intersection disappears, since the orbital plane becomes 
the invariable plane and the associated inclination an- 
gle vanishes. In this case, the related angles T and <p 
are not individually defined, though Eq. (4.27) gives us 
<p + T = 0, as required. We emphasize that the related 
differential equations, given by Eqs. (4.29) and (4.32), are 
obtained using 6^0, corresponding to S ^ 0. There- 
fore, Eqs. (4.29) and (4.32) are not defined if = 0. 
Further, the apparent non-vanishing of the time evolu- 
tion of T in the non-spinning limit is also attributable 
to the similar behavior for iv p . However, a close scrutiny 
of Eqs. (4.19), (4.22) with (4.23), (4.27) and (4.28) for 
5^0, which implies 0^0 and 6 — > tt/2, reveals that 
the binary dynamics in the non-spinning case is describ- 
able in terms of r and <fi, where <fi = ip + T , as explained 
earlier. 

Our parametric solution generalizes a restricted anal- 
ysis considered in Ref. [28]. The analysis, given in 
Ref. [28], neglects the precessional motion of the spin 
vectors and restricts S e s, denoted by C m Ref- [28], to 
lie along ez. Due to these restrictions, that analysis 
provides only a parametric solution to r. We have, in 
the other hand, consistently taken into consideration all 
the leading order spin-orbit interactions and obtained a 
parametric solution to the entire binary dynamics. In 
a future communication we will connect the above pre- 
sented parametric solution to quantities related to obser- 
vation, especially associated with binary pulsars [45]. 

In the next section, we will combine the Keplerian- 
type parametrization developed above with the one pre- 
sented in Sec. III. 



V. THIRD POST NEWTONIAN ACCURATE 
GENERALIZED QUASI-KEPLERIAN 
PARAMETRIZATION FOR COMPACT 
BINARIES IN ECCENTRIC ORBITS WITH THE 
FIRST-ORDER SPIN-ORBIT INTERACTION 



We are now in a position to have parametric solution 
to the dynamics, defined by the Hamiltonian H, as given 
by Eq. (2.2). This is done by combining consistently 
the parametrizations presented in the previous two sec- 
tions, Sec. Ill and Sec. IV. The two parametric solutions 
of Sec. Ill and Sec. IV can be added linearly to obtain 
the PN accurate binary dynamics, given by Eqs. (2.2) 
and (2.3), for the following reasons. Recall that in the 
Hamiltonian, given by Eqs. (2.2) and (2.3), the spin- 
orbit contribution is added linearly to the PN accurate 
non-spinning contributions. Moreover, as we are con- 
sidering only leading order spin-orbit interactions, spin- 
orbit contributions only cross with the Newtonian terms 
in the Hamiltonian. This is why we can treat the spin- 
orbit contributions and the non-spinning PN contribu- 
tions separately and later add the two parametric solu- 
tions linearly. However, as a cautionary note, we state 
that care should be taken while merging the above two 
parametric solutions to avoid adding similar contribu- 
tions twice. Finally, we display below, in its entirety, 
the parametric solution to the conservative third post- 
Newtonian dynamics of spinning compact binaries mov- 
ing in an eccentric orbit, in ADM-type coordinates, with 
spin effects restricted to the leading order spin-orbit in- 
teractions: 



r = r cos ipi + r sin tp j 

L = Lk, 

S = Jez — Lk , 



(5.1) 
(5.2) 
(5.3) 



where the basic vectors (i,j,k) are explicitly given by 

i = cos Te.Y + sin Tey , (5.4a) 

j = — cos sin Tex + cos cos Tey + sin ez , 

(5.4b) 

k = sin sin Tex — sin cos Tey + cos ez , (5.4c) 



and 



r = a r (1 — e r cos u) , 

t /i i \ • . (9u . 56* \ , v . / fit . fet 

I = n [t — to) = u — e t sin u + y— j- + — g- J (v — u) + I + -^r 



— 7t sin 2v H 7t- sin 3v , 

c 6 c 6 



(5.5a) 
(5.5b) 



if - = (1 + k)v + I i£ + ^ J sin 2v + (£J + ^) sin 3v + ^ sin Av + ^ sin 5v , (5.5c) 
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T - T = ^2°r3 ( v + es ^ nv ) ■ 
C Li 

where v = 2 arctan 



1 - e a 



1/2 



tan ■ 



(5.5d) 
(5.5e) 



The post-Newtonian accurate expressions for the orbital elements e 2 , k, and e 2 and the post-Newtonian 

orbital functions g± t , get, fit, ht, Ht, h & t, fi V , fe v , 9i v , 36 V , «6 V , and h 6ip , in terms of E, L, S, ij and a, read 



94t 
96t 

fit 
fet 



(-2E) 



(-2E) 



192c 6 



4c 2 

3 — 977— 6^ 2 + 3?7 3 + 



7 + 7? + 4x so cos a- 



16c 4 



1 + 10r? + rf + 



1 



{-2EL 2 



■(-68 + 44??) 



1 



{-2EL 2 ) 



(864 - 3tt 2 7? - 2212r/ + 432?? 2 



2 Z L2 ^ 2 (-6432 + 1348877 - 240tt 2 t7 - 768?? 2 ) J 



(5.6a) 



= 1 + 2EL Z 



{-2Ef 



(-225) J 
Ac 2 \ 



24 - 477 - 5(3 - ri)(-2EL 2 ) - 16(1 + EL 2 ) Xso cos a 



8c 4 
(-2E) 3 
192c 6 
4 

{-2EL 2 ) 



52 + 277 + 2?? 2 - (80 - 55r? + 4?? 2 ) (-2EL 2 ) + 



L 

(17-IItj) 



(-2EL 2 ) 

-768 - 6tt 2 t? - 344?? - 2I677 2 + 3{-2EL 2 ) (-1488 + 1556?? - 319r? 2 + 4?? 3 

192 



- 821277 + 177tt 2 77 + 48O77 2 ' 



(-225) 



(-2Ef' 2 ^l+ x - 



{-2Ef 



1 

3072c 6 
5760 

V-2EL 2 

1 + 2EL 2 + 



8c 2 v " 128c 4 

-29385 - 4995?7 - 31577 2 + 13577 3 

(17 — 977 + 2t7 2 

-2E 
~Ac 2 ~ 



{-2EL 2 ) 2 
555 + 30??+ II?? 2 
16 



(134 - 281?? + 5tt 2 77 + I677 2 ) | , 
(5 - 27?) 



(5.6b) 



192 



{-2EL 2 f/ 2 



8 + 877 + (17 - 7r])(-2EL z ) - 8x so cosa- 



V-2EL 2 
(10080 + 123tt 2 77 - 139527? + 1440?? 2 



S 



(5.6c) 



^ 2E ) (8 + 477 + 2O77 2 - (-2252, 2 )(112 - 47?7 + I677 2 ) + 2A\/-2EL 2 (5 - 277) 
8c 4 ( 



{-2EL 2 ) 
{-2Ef 



(17-11??) 



24 



V-2EL 2 



(5 - 27?) 



192c 6 



1 24 (-2 + 5?;) (-23 + 10t? + 4?? 2 ) - 15 (-528 + 200?/ - 77?? 2 + 24?? 3 ) (-2EL 2 ) 



72(265 - 193t? + 46t? 2 ) -2EL 2 



{-2EL 2 ) 



(6732 + 1177T 2 ?? - 12508?? + 2004?? 2 



===== (16380 - 19964?? + 123tt 2 t? + 3240?? 2 ) - - 2 EL 2 ) 3 / 2 ( 10080 + 123?r2? 7 _ 13952?? + 1440?? 2 ) 



+ { _^L 2 ) 2 (134 " 281V + ^ + 16??2) } 
3(-2£) 2 5 - 2?? 



V-2EL 2 ' 



{-2Ef 



1 



192 \ (-2EL 2 ) 3 / 2 

_1 ( -225) » A / 1+ 9,F?r.2 

8 V-2EL 2K 
(-2E) 3 1 



(10080 + 123tt 2 7? - 13952?? + 1440?? 2 



(5.6d) 
(5.6e) 

-==L= (-3420 + 1980?? - 648?? 2 ) | , (5.6f) 

(5-6g) 



192 VI + 2EE 2 I (-2EL 2 ) 3 / 2 



■(1728 - 414877 + 3tt 2 7? + 600?? 2 + 33?? 3 ) + 3v /Z 2£T 2 ??(-64 - 4t? + 23?? 2 ) 
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1 



«6t 
k 



96<p 

4 



y/-2EL 2 
(-2E) 3 (1 + 2EL 2 



32 (-2EL 2 ) 3 / 



(-1728 + 4232?7 - 3tt 2 7? - 627t/ 2 - 1057y 3 ) 
77(23 + 12? ? + 6t/ 2 ) , 



13(-2£) 3 ^ 3 



192 



1 + 2EL< 
-2EL 2 



2\ 3/2 



<z 2 L 2 



-g Xso cosa-^ 



S , ( -2g) 
4c 2 



-5 + 2?] + 



35-lOrA (-2E) 2 
-2EL 2 J + 384c 4 



120 - 120?7 + 96t7 2 



1 



-2EL 2 ) 
1 



{-2EL 2 ) 2 
(-2E) 2 (1 + 2EL 2 ) 



(-10080 + 13952t/ - 123tt 2 t7 - 1440t/ 2 ) 
(36960 - 40000?? + 615tt 2 ?7 + I68O77 2 
77(1 - 3 ?7 ) , 



(-2EL 2 ) 2 



1 (-2£) 3 f 



256 (-2EL 2 ) \ 
1 



(44 + I6O77 — 96?7 2 )?7 



(-2EL 2 ) 2 
3(-2£) 2 



{-2EL 2 ) 

(256 + 497r 2 77 - 98O77 - 672? ? 2 - 4O77 3 ) 



-256 - 49tt 2 77 + 109677 + 624?/ 2 - 80t/ 3 



>r 



32 (-2EL 2 ) 2 



(1 + 2EL Z ) 



2\3/2 



(-2E) 3 VTT2EL 2 
768 (-2EL 2 
1 



77^ -27t7 + 78t7 2 



1 



(-2EL 2 



-(220 + 3tt 2 + 967/+ 1507/ 2 ) 



{-2EL 2 ) 2 
(-2E) 3 (1 + 2EL 2 ) 2 



128 {-2EL 2 ) 3 
h(-2E) 3 rj 3 
256 (-2EL 2 ) 3 

2 , i-m 



(220 + 3tt 2 - 1207/ + 45?7 2 ) 
77(5 + 2877 + IO77 2 ) 



1 + 2EL 

(-2£) 2 

16c 4 
(-2B) 3 



(1 + 2EL 2 ) 5/2 , 
24 - (15 - r/)(-2EL 2 ) + 8(1 + 2EL 2 ) Xs 



4c 2 



- 32 + 1767? + I877 2 - (-2t3L 2 )(160 - 30?7 + 3?7 2 



8(3 + 4£L 2 )x so cosa- 



(408 - 232?; - 15?/ 2 



(-2EL 2 ] 



+ 



l - 16032 + 276477 + 3tt 2 t) + 4536t7 2 + 234?/ 3 - 36 (248 - 80?/ + 13?; 2 + 77 s ) (-2EL 2 ) 
(2456 - 26860?/ + 581tt 2 ?7 + 2689?/ 2 + IO77 3 ) 
n 2-2 (27776 - 65436?7 + 13257T 2 ?/ + 3440?/ 2 - 70?/ 3 ) 

LsIj Li ) 



384c 6 
6 

[-2EL 2 ) 
3 



(5.6h) 
(5.6i) 

(5-6j) 



(5.6k) 
(5.61) 

(5.6m) 
(5.6n) 

(5.6o) 
(5.6p) 
(5.6q) 



(5.6r) 



We recall that e, appearing in Eq. (5.5d), is given by the Newtonian contribution to e r . The three eccentricities e r , 
et and e v are related to each other by post-Newtonian corrections 
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(5.7b) 



and hence it is possible to describe the binary dynamics 
in terms of one of the eccentricities. 

We emphasize that, while adapting the parametriza- 
tion to describe the binary dynamics, care should be 
taken to restrict orbital elements to the required post- 
Newtonian order. 

As noted earlier, 6, the precessional angle for k, is 
given by 



sin 6 



cos 6 



S'sina 

L + S cos a 
J 



(5.8a) 
(5.8b) 



where the magnitude of the total angular momentum is 
given by J = (L 2 + S 2 + 2LScosa) 1 / 2 . 

Finally, we want to emphasize that the parametric so- 
lution, given by Eqs. (5.1)-(5.8), describes the entire con- 
servative post-Newtonian accurate dynamics of a spin- 
ning compact binary in an eccentric orbit, when the spin 
effects are restricted to the leading order spin-orbit in- 
teraction. This means the parametrization consistently 
describes not only the precessional motion of the orbit in- 
side the orbital plane, but also the precessional motions 
of the orbital plane and the spins themselves. 

It is possible, though tedious, to show that the com- 
bined parametric solution, given by Eqs. (5.1)-(5.8), is 
indeed a parametric solution to the total orbital dynam- 
ics prescribed by Eqs. (2.2) and (2.3). To show that, we 
have computed, using the combined parametric solution, 
PN accurate expressions for r, 9 and <j> in terms of E, 
L, L ■ S eff and r with the help of Eqs. (4.27) and (4.28). 
These expressions were found to be in total agreement 
with PN accurate expressions for f, 9 and <f> computed 
from the total (reduced) Hamiltonian, given by Eqs. (2.2) 
and (2.3), using the Hamilton's equations of motion [See 
Eqs. (4.7)-(4.11) and the related discussions]. 



VI. THE LEADING PART OF THE 
QUADRUPOLAR GRAVITATIONAL-WAVE 
POLARIZATIONS FOR SPINNING COMPACT 
BINARIES IN ECCENTRIC ORBITS 

As an application of the above parametrization, we ob- 
tain, for the first time, explicit expressions for h + and h x , 
suitable to describe gravitational radiation from spinning 
compact binaries moving in eccentric orbits, correspond- 
ing to the cases considered in Sec. IV. The gravitational 



wave polarization states, h + and h x , are usually given 

by 



\ (PiPj - QiQj) h 
1 , 

- {piqj +Pjqi)h 



TT 

ij 



TT 

ij 



(6.1a) 
(6.1b) 



where p and q are two orthogonal unit vectors in the 
plane of the sky, i.e. a plane transverse to the radial 
direction AT (the line of sight), linking the observer to 
the source (see Fig. 2). 

The transverse-traceless (TT) part of the radiation 
field, hJ 3 T , which depends on the dynamics of the com- 
pact binary, is expressible in terms of a post-Newtonian 
expansion in (v/c). Symbolically, the TT radiation field 
may be written as 



lTT 



1 



(0) 



„5 "»j 



1 



(2) 



1 



(3) 



(6.2) 



In this paper, for simplicity, we will restrict hj^ to its 
leading 'quadrupolar' order, namely to /c 4 , and de- 
note it by hJ?\Q. However, higher PN corrections to 
are available in the literature. The 2PN corrections 

l j 

to hj^ for compact binaries moving in general orbits are 
given in Refs. [25, 46]. The spin-orbit and spin-spin in- 
teractions directly contribute to hJ 3 T and were computed 
in Ref. [23]. For compact binaries, in circular orbits, 
the explicit expressions for 2.5PN accurate hT 



TT 



cor- 



rections that include /c 9 in harmonic coordinates — 
were recently derived in Ref. [5]. We recall that post- 
Newtonian corrections to hj? are usually given in har- 
monic coordinates, which differ from the one we employed 
here at 2PN and 3PN orders. However, using coordinate 
transformations that link harmonic and ADM-type co- 
ordinates, given in Ref. [12], it is possible to obtain post- 
Newtonian corrections to the amplitudes of h+ and h x 
in ADM-type coordinates. 

The explicit expression for hJ } T , in the leading part of 
the 'quadrupolar' approximation, reads 



lTTI 

"bi 



ijkrn 



(AT) 



GM 



(6.3) 



where Vijkm(N) is the usual transverse traceless pro- 
jection operator projecting vectors normal to AT, where 
N = R'/R' is the line-of-sight unit vector from the 
observer to the binary, and R' — \R'\ is the corre- 
sponding radial distance. We also used := ViVj, and 
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Figure 2: The convention we adopted to link the orbital frame 
(i,j,k), the invariable frame (ex,ey,ez) and the frame 
(p,q,N) associated with the observer. Since we choose the 
line-of-sight unit vector N to lie in the ey-e^-plane, we may 
align the polarization vector p along ex, where the plane of 
the sky meets the invariable plane. This implies that q also 
lies in the ey-ez-plane. Therefore the frames (p,q,N) and 
(ex , ey , ez) are connected by a constant angle ia, which is 
the constant inclination angle between N and J. The po- 
larization vectors p and q span the plane of the sky. The 
inclination of this plane with respect to the orbital plane is 
the orbital inclination i. The inclination of the orbital plane 
with respect to the invariable plane is denoted by the constant 
angle O. 



riij := n.irij, where Vi and rii are the components of the 
velocity vector v = dr jdt and the unit relative separation 
vector n = r/r, where r = |r|. 

Using Eqs. (6.1) and (6.3), we obtain /i+|q and /i x |q, 
the expressions for gravitational wave polarizations when 
their amplitudes are restricted to the leading quadrupolar 
order, for compact binaries moving in general orbits as 



h+\ 



2Gfi \, , / GM 

-T^j \ iPiPi - qiOj) I «« - — i 

-^[(p-n) 2 -( g .n) 2 ]|, 



(6.4a) 



2G[i 
4Gfi 



GM 



GM 

(j> ■ v)(q ■ v) (p ■ n)(q ■ n) 



(6.4b) 

The expressions for /i+|q and h x |q, for spinning com- 
pact binaries in eccentric orbits are obtained by adapt- 
ing the radial separation and velocity vectors, r and 
v = dr/dt, to the orthonormal triad (p,q,N). This is 
easily achieved in the following way. We deduce that 
the explicit parametric representation for r, given by 
Eqs. (5.1) and (5.4), in terms of the reference triad 
(ex, ey, ez) is given by 

r = r(cos T cos <p — cos O sin T sin ip)ex 
+ r (sin T cos ip + cos 9 cos T sin ip)ey 
+ r sin sin ip ez ■ (6-5) 

In the above equation, the angles ip and T along with 
r depend on the orbital dynamics. However, 6 is a 
constant angle — a direct consequence of choosing two 
distinct cases for the parametrization (see discussion in 
Sec. IV). 

To compute /i+|q and /i x |q, one needs to choose a 
convention for the direction and orientation of the orbit 
with respect to the plane of sky. In our approach, the 
orthonormal triad (p, q, TV) is connected to the reference 
triad (ex,ey,ez) by a constant angle «o (see Fig. 2). 

In particular, we connect the reference triad 
(ex,ey,ez) to the observer triad (p,q,N) by a rota- 
tion around p via a constant angle Iq 



(6.6) 



where Ci and Si are shorthand notations for cos iq and 
sinio. It is now straightforward to get the explicit ex- 
pressions for the radial separation and velocity vectors in 
terms of the orbital triad (p,q,N), associated with the 
observer. The vectorial equation for r reads 

r = r(cos T cos ip — Cq sin T sin tp)p 

+ r [d sin T cos <p — (Si Sq — Ci Cq cos T) sin <p] q 

+ r [Si sin T cos ip + (Ci Sq + Si Cq cos T) sin ip] N , 

(6.7) 

where Cq and Sq are shorthand notations for cos 9 and 
sin 6. The radial velocity vector v = dr/dt is given by 
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cos ip 
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+ | — rCi Q (ip + CeT) sin T + Cj Cef cos T — Si S@r simp 
+ Cj fsinT + rCi (C®ip + T) cosT — Si S®r(p cos^jqr 
+ | —rSi ((p + CeT) sinY + S^Cef cos T + C ig Sef simp 
Si r sin T + rSi {C®<p + T) cos T + Ci S®rip cos AT 



(6.f 



Note that while computing v, we have kept and io as constant angles. We point out that the way we defined the 
angle io is similar to the manner in which the 'inclination' was defined in Ref. [47] and Ref. [23] [see, e.g., discussions 
prior to Eqs. (4.23) in Subsec. C under Sec. IV in Ref. [23]]. 

It is now straightforward, though lengthy, to obtain the explicit expressions for gravitational wave polarizations for 
spinning compact binaries moving in eccentric orbits, when the dynamics includes the leading spin-orbit interaction. 
The expressions, H + \q and h x |q, giving the gravitational wave polarizations with amplitudes due to the quadrupolar 
contributions, read 
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[GM 
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C io S'|rrtcos2T 



(6.9b) 



We would like to remind (again) the reader that the or- 
bital phase is denoted by ip, the angle T describes the 
phase of the line of nodes i, (p = dtp/dt, Y = dT/dt, 
and r = (n-v). The parametric PN accurate expressions 
for these time derivatives can easily be computed using 
the parametric solution as r = dr/dt = (dr/du)(du/dt), 
p = dip/dt = (dip/dv)(dv/du)(du/dt) and T = dT/dt = 
[dT /dv){dv/du){du/dt). 

Notice that the temporally evolving (observational) or- 
bital inclination i (see Fig. 2), defined by cosi = N ■ k, 
does not enter the expressions for /i+|q and /i x |q. How- 
ever, using io, O and T, the time evolution of i is simply 
given by 



cos i = TV ■ k = d Cq 
sini = 17V x k\ 



Si S@ cos T . 



(6.10a) 



[(S l0 C e + C io S@ cosT) 2 + S| sin 2 T] 
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(6.10b) 



The equation for cosi is in agreement with Eq. (6) in 
Ref. [47]. We may also obtain a differential equation for 
the time evolution of i, as given in Ref. [35]. To get that, 
we differentiate cost = N ■ k with respect to time and 
making use of k = L/L with Eq. (4.2a) and Eq. (6.10b). 
In this way we deduce the following differential equation 
for the orbital inclination i 



di 



1 TV x k 

rSel ' — ; — — 



1 TV x k 

eff ' 17V x fcl 



(6.11) 



which agrees with Eq. (5.15) of Ref. [35]. 

The Eqs. (6.9) are useful for scenarios where the angle 
is a constant. This is true for the two cases where 
our parametric solution is applicable (see discussions in 
Sec. IV after Eqs. (4.22)). 

The approach to compute H+\q and H x \q may be 
adapted to include higher order spin effects like spin-spin 
interactions. In this case, 6 will no longer be a constant 
angle and expressions for the velocity vector v and hence 
/i+|q and H x \q will also depend on O. 

The non-spinning limit, as given by Eqs. (6) in 
Ref. [24], is obtained in the following way. We note that 
when S — > 0, J is identical to L, which implies 6^0, 
9 — > it/ 2 and io — ► i, i being the angle between TV and 
L. In this limit, Eqs. (4.27) and (4.28b) indicate that 
T + p = 4> and T + ip = <fi. This leads to Eqs. (6) in 
Ref. [24] which reads 
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where C, and Si are shorthand notations for cosz and 
sini. 

The temporal evolution of h + and h x , given by 
Eqs. (6.9), may be obtained by adapting the 'phasing 
formalism', presented in Ref. [24], and this will be re- 
ported soon [48]. However, it should be noted that the 
parametrization, given by Eqs. (5.1)-(5.8), will be suf- 
ficient to determine the 3PN accurate conservative time 
evolution, which includes the leading order spin-orbit in- 
teraction. 



VII. CONCLUSIONS 

We have presented Keplerian-type parametrization for 
the solution of post-Newtonian accurate conservative dy- 
namics of spinning compact binaries moving in eccentric 
orbits. The above PN accurate dynamics consisted of 
3PN accurate conservative orbital dynamics, associated 
with non-spinning compact objects, influenced by the 
leading order spin-orbit interactions. The orbital ele- 
ments of the representation were explicitly obtained in 
terms of the conserved orbital energy, angular momen- 
tum and a quantity that characterizes the leading order 
spin-orbit interactions in ADM-type coordinates. Our 
parametric solution is applicable in the following two dis- 
tinct cases, namely, case (i), the equal-mass-case and 
case (ii), the single-spin-case. We also derived expres- 
sions for gravitational wave polarizations suitable to de- 
scribe gravitational radiation from spinning compact bi- 
naries moving in eccentric orbits. 

The present work has several possible applications and 
some of them are currently under investigation. Em- 
ploying the 'phasing formalism', given by Ref. [24], along 
with our parametric solutions and expressions for h + and 
/i x , we will, in the near future, obtain gravitational wave 
polarizations suitable to describe gravitational radiation 
from spinning compact binaries moving in inspiralling ec- 
centric orbits [48]. Our parametric solution for the two 
cases considered implies that the associated PN accurate 
conservative binary dynamics will not be chaotic. This 
is so as our solution can analytically determine the asso- 
ciated dynamics. It will be interesting to explore (again) 
numerically the scenarios, where our parametrization is 
valid, and investigate if the numerical solutions still pre- 
dict chaos. Finally, as mentioned earlier, the fully 2PN 
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accurate 'timing formula' may be derived using our para- 
metric solution as a crucial input. This may be required 
as spin-orbit interactions essentially enter at 2PN order 
for binary pulsars and their effects leave observational 
signature [35, 49-51]. 

It is highly desirable to extend the present work in the 
following directions. In the literature, there exists a PN 
accurate orbital equations of motion for spinning com- 
pact binaries, where spin-orbit interactions are also PN 
accurate [22]. It will be interesting to see if a parametric 
solution is possible for the above dynamics. Naturally, 
it will be interesting to include spin-spin effects in to 
our parametric solution. However, all these extensions 
are not straightforward and the parametric solutions for 



these cases, most probably, will not be as elegant as the 
one presented in this paper. 
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